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Introduction 

In  the  following  sections  we  discuss  the  scientific  accomplishments 
made  daring  the  course  of  the  19  1/2  months  of  this  contraction  the 
work  was  in  the  area  of  gasdynamics,  and  our  main  area  of  research 
has  been  cases  of  flows  containing  shockwaves.  For  the  most  part  we 
have  taken  an  approach  which  straddled  the  borderline  of  analytical 
and  numerical  methods.  This  is  best  exemplified  by  the  work  discussec^: 
in  Sec.  l<^in  which  classical  boundary  layer  methods  are  used  to  develop 


variable  mesh  numerical  schemes^  and.  in  Secs.  2 and  3*in  which  the 
theory  of  characteristics,  Riemann  invariants,  and  shock  expansion 
theory  are  all  used  to  obtain  a robust  method  of  numerical  integration 
which  is  orders  of  magnitude  faster  than  any  other  known  scheme.^ 

Daring  the  contract  period  the  four  following  manuscripts  haVe 
been  prepared: 

1.  "A  Variable  Mesh  Finite  Difference  Method  for  Solving  a Class  of 
Parabolic  Differential  Equations." 

"An  Approximate  Solution  in  Gasdynamics . " 

"On  the  numerical  Integration  of  the  Gasdynamic  Equations." 

"On  the  Equations  Governing  a Trailing  Wake." 

Of  these  the  first  has  been  published  in  the  SIAM  Journal  on 
Numerical  Analysis,  and  the  other  three  have  been  submitted  for  publication. 
Two  additional  manuscripts  are  in  preparation: 

1.  "Two-Dimensional  Supersonic  Flow:  I.  Analytical  Approximation." 

2.  "Two-Dimensional  Supersonic  Flow:  II.  Exact  Numerical  Integration." 


2. 

3. 

4. 
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In  addition  to  the  work  covered  by  these  manuscripts,  several  other 
projects  were  undertaken  but  were  not  felt  to  warrant  a publication. 

These  are  also  discussed  below. 
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Numerical  Integration  of  Equations  of  the  Navier-Stokes  Type 
A number  of  situations  in  fluid  mechanics  lead  to  flows  containing 


"boundary  layers"  where  flow  variables  experience  rapid  variation. 

In  the  case  of  shock  waves  a priori  location  of  this  region  is  not  known 
(this  is  also  true  of  the  wake  roll-up  region  behind  a moving  body) . 

The  numerical  integration  of  the  fluid  equations  in  such  situations 
present  a number  of  difficulties.  In  any  real  situation,  where  Reynolds 
numbers  are  large,  transition  layers  are  small  and  a numerical  scheme 
based  on  a uniform  mesh  is  unfeasible.  To  deal  with  this  situation 
we  have  investigated  a variable  mesh  scheme  for  the  equation 


(1) 


rLu  m 


F(X,  t,  U,  U , u ) 


From  the  mathematical  viewpoint  this  represents  a parabolic  partial 
differential  equation  and  it  is  known  that  such  equations,  even  for 
smooth  data,  give  rise  to  region6  .of  rapid  variation  for  which  a priori 
error  bounds  fail. 

A particular  example  of  equation  (1)  is 


(2) 


kt-n/2u|u  |fu 
bt  ox  2 ^2 


which  for  n «-l  is  the  -equation-of  Burgers  (1) . ■'More  recently  Eqn.  (2), 
which  we  call  the  modified  Burgers  equation,  has  been  derived  from  the 
Navier-Stokes  equations  for  problems  involving  explosions  and  sonic 
booms  [2],  (3}. 

He  have  developed  a method  for  dealing  with  problems  of  this  type. 


4 


which  is  uniformly  second  order  accurate.  A variable  mesh  (in  x)  is 
generated  based  on  boundary  layer  methods.  The  neighborhood  of  the  shock 
contains  a zone  of  high  density  grid  points,  and  a low  density  grid 
is  used  in  regions  of  slow  variation.  The  grid  system  is  rezoned  during 
the  course  of  the  numerical  integration,  i.e.,  grid  generation  is  an 
integral  part  of  the  numerical  code.  Further  reduction  of  the  numerical 
problem  is  obtained  from  use  of  approximate  shock  trajectories,  which  can 
be  derived  from  asymptotic  theory. 

The  results  of  this  investigation  extend  by  three  orders  of  magnitude 

a similar  calculation  (using  a uniform  mesh)  found  in  [4] . As  an  example 

r -5 

if  n - 2 and  the  reciprocal  Reynolds  number  * is  10  , a uniform  mesh 

approach  requires  > 105  mesh  points.  For  the  method  discussed  here 
600  points  are  required  at  the  initial  instant.  With  the  course  of  time 
this  decreases  steadily,  e.g.,  at  t - 2000  only  300  points  are  required. 

The  following  figures  give  in  graphical  form  the  evolution  of  an  initial 
discontinuity  for  the  three  Reynolds  numbers,  R ■ $ * ■*  10^,  10*,  10^. 
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A report  of  this  work  appears  in  the  following  publication: 

T.1I.  Chong,  A Variable  Mesh  Finite  Difference  Method  for  Solving 
a Class  of  Parabolic  Differential  Equations  in  One  Space  Variable,  SIAM  J. 
Numer.  Anal.  15^  4 (August  1978) ,pp.  835-857. 

Abstract.  A variable  mesh  finite  difference  scheme  for  a class  of 
parabolic  differential  equations  which  exhibit  shock-like  structures  is 
developed.  It  is  shown  that  a properly  chosen  variable  mesh  will  yield 
results  comparable  in  accuracy  to  one  using  a much  finer  uniform  mesh. 
Computable  criteria  and  schemes  for  generating  such  variable  meshes  are 
given.  A scheme  is  then  applied  to  the  Burgers'  and  modified  Burgers' 
equations  with  a small  viscosity.  Excellent  agreement  is  obtained  with 
known  exact  solutions. 

Further  work  should  consider  multiple  shocks,  shock  interaction,  and 
higher  dimensional  discontinuities. 

References 

[1]  J.M.  Burgers,  A mathematical  model  illustrating  the  theory  of 
turbulence,  Advances  in  Appl.  Mech.,  1 (1948),  pp.  171-199. 

12)  T.H.  Chong  and  L.  Sirovich,  Nonlinear  effects  in  supersonic 
dissipative  gasdynamcis.  Part  2,  J.  Fluid  Mech.,  58  (1973),  pp.  53-63. 

13]  L.  Halabisky  and  L.  Sirovich,  Evolution  of  finite  disturbances  in 
dissipative  gasdynamics,  Phys.  Fluids,  16  (1973),  pp.  360-368. 

(4)  P.L.  Sachdev  and  R.  Seebass,  Propagation  of  spherical  and  cylindrical 
N-waves,  J.  Fluid  Mech.,  58  (1973),  pp.  192-205. 

(5)  L.  Sirovich  and  T.H.  Chong,  Supersonic  flight  in  a stratified 
sheared  atmosphere,  Phys.  Fluids,  17  (1974),  pp.  310-320. 


2,  Unsteady  Gasdynaroics — Analytic  Approximation  and  Numerical  Integration 


We  have  pursued  a novel  method  of  treating  gasdynamic  flows  in 
which  shock  waves  occur.  Our  research  has  produced  a new  analytical 
approximation  and  a new  numerical  procedure  for  the  exact  integration 
of  gasdynamic  flows  with  shocks.  As  a result  of  the  analytic  approach 
we  have  been  able  to  recast  the  gasdynamic  equations  in  a form 
which  we  believe  is  optimal  for  machine  integration.  The  result  has  been 
a highly  accurate  and  extremely  rapid  numerical  pro  -dure.  Unlike 
existing  methods,  where  shock  waves  oscillate,  and  often  fail  to  stabilize, 
the  present  method  is  able  to  fix  a shock  wave  once  and  for  all. 

We  consider  the  situation  as  depicted  in  the  figure. 


A piston  is  moved  into  a gas  at  rest  with  the  trajectory  P (t) . Leading 
and  trailing  shocks  X^(t)  and  X^(t)  are  formed.  In  brief  the  analysis 
is  as  follows: 

Riemann  invariants  (normalized) 


■ , * 2c 

(1)  r.Oj^ 


are  introduced.  ((1)  is  for  an  ideal  gas,  the  modification  of  the  develop- 
ment for  non-ideal  gases  is  relatively  minor.)  A new  coordina±e  sys item 
( , p ) is  introduced  through 

a*  **  A 

— + u r—  = 0 
dt  dx 


SI  * (o  ♦ c>  ■ 0 


The  gasdynamic  equations  then  take  on  the  form 


dS  A 

if ' 0 


<3>  &r  ‘ s sf 

Ml  Ff'v.'1* 

The  time,  t,  is  governed  by 


Vj  + ln  c,t*  " 


and  the  spatial  variable,  x,  can  then  be  determined  from 


V.  " <u  + c^t0. 


• ' 4 • - re-  ••  • 

• ’ ' J»  ■ %r*  *• ' _ N 


These  equations  are  to  be  augmented  by  shock  and  boundary  conditions, 
which  jure  not  repeated  here. 

The  first  step  in  the  procedure  follows  from  an  approximation 
which  occurs  in  shock  expansion  theory  [1] , [2),  [3],  viz., 

(6)  r r (<* ) 

o 

With  adoption  of  (6)  we  drop  equation  (4).  Eqs.  (2),  (3),  (5),  and 
(6)  constitute  the  approximate  theory  and  cam  be  solved  in  terms  of 
quadratures  and  the  functions  S(tf-),  rQ(rt  ) which  are  evaluated  at  the 
shock  wave.  The  results  of  this  approximation  are  remarkably  accurate. 
Even  in  the  neighborhood  of  ionizing  shocks  the  approximate  and  exact 
solutions  differ  by  less  than  one  percent. 

The  exact  numerical  procedure  uses  the  approximate  solution  as  a 
first 'step  in  an  iterative  procedure.  The  approximate  solution  is 
substituted  into  the  right  hand  side  of  (4)  and  a new  value  of  r , IS  ) 
obtained.  This  in  turn  is  used  to  solve  (2),  (3),  (5),  and  (6),  and  the 
iteration  repeated  if  necessary.  According  to  the  present  version 
of  the  numerical  codes,  this  procedure  converges  in  less  than  three 
iterations. 

As  an  example  which  is  typical  of  the  integration  procedure, 
we  consider  a parabolic  piston  motion, 

P(t)  - MQtU  - t/2) 

where  M is  the  Mach  number  of  the  piston  at  the  initial  instant.  At 


* 3 the  flow  is  considerably  away  from  equilibrium.  For  example, 
the  temperature  behind  the  initial  shock  is  four  times  its  initial  value. 
The  result  of  the  numerical  integration  for  MQ  ■ 3 is  shown  in  the 
following  four  figures. 

The  first  of  these  figures  indicates  the  particle  paths,  ^ equal 
+ * 

constant,  the  c characteristics,  p equal  constant,  and  the  shock  waves, 
heavy  lines,  in  the  physical  plane.  Each  of  the  following  three  figures 
gives  the  variation  of  the  physical  quantities,  S,  r , p,  c on  either 
particle  paths  or  C+  characteristics.  Thus  the  solution  to  the  piston 
problem  is  furnished  entirely  in  graphical  form. 
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The  results  describe  in  this  section  are  reported  on,  in  detail, 
in  the  following  two  manuscripts: 

L.  Sirovich  and.T.H.  Chong,  An  Approximate  Solution  in  Gasdynamics . 
T.H.  Chong  and  L. ' Sirovich,  On  the  Numerical  Integration  of  the 
Gasdynamic  Equations. 

Both  of  these  manuscripts  have  been  submitted  for  publication. 

•'  References 

tl)  W.D.  Hayes  and  R.  Probstein,  Hypersonic  Flow  Theory  (2nd.  Ed.), 
Academic  Press,  New  York,  (1966),  p.  497ff. 

(2)  J.J.  Mahony,  J.  Aeronat.  Aci.  22,  673-680  (1955). 

R.  Meyer,  Quart.  Appl.  Math.  1£,  433-436  (1957). 
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3.  Steady  Supersonic  Flows 

In  parallel  with  the  treatment  described  in  Sec.  2,  we  have  also 
investigated  the  case  of  two-dimensional  supersonic  flows.  The  analysis 
and  methods  summarized  below  applied  to  flows  ranging  from  the  low 
supersonic  limit  to  high  supersonic  speeds. 

The  two-dimensional  form  of  the  gasdynamic  equations  in  characteristic 
form  follows  from  the  introduction  of  the  Prandtl  angle 


/<.  “ sin 


and  the  flow  angle 

0 * tan  1 £ 

as  dependent  variables.  Then  with  the  introduction  of  the  Riemann 
invariants 


(1) 


r ± « 0 + p(u) 


where  the  Prandtl  function  P is  defined  by 


the  governing  equations  are: 


(3)  ||  - 0 


<«>  ri  • $ *in  A4  001 V d* 


A 


% 


In  these  the  physical  plane  (x,  y)  has  been  transformed  to  (#<,/)- 
plane  by  means  of 

(6)  y^  “ tan  0 x^ 

(7)  y^  « tan  ( 0 +yw.  )x^ 

+ 

The  first  of  these  corresponds  to  streamlines  and  the  second  to  C characteristics. 

The  method  of  solution  (in  the  upper  half  plane)  follows  from 
what  we  refer  to  as  the  approximat  ion. of  shock  expansion  theory  (1-3) , 
viz.  that  r is  "almost  constant"  on  streamlines.  We  therefore  replace 
Eq.  (6)  with  the  condition 

(8)  r » r (rt  ) 

Both  r (/)  and  S(e<  ) are  to  be  determined  at* the  -shock. 

The  solution  to  the  approximate  theory  as  given  by  (3),  (4),  (8), 
as  well  as  shock  relations  not  listed  here  leads  to  solution  in  terms  of 
quadratures.  Furthermore,  it  leads  to  results  at  the  body  which  are 
within  one  percent  of  the  exact  results  (to  be  described  below) , even 
at  the  highest  Mach  numbers. 

An  exact  numerical  integration  also  proceeds  from  the  above  formulation. 

The  results  of  the  approximate  theory  are  substituted  into  the  right  hand 
side  of  Eq.  (5)  which  leads  to  a new  determination  of  r (#(  ) (note 

that  r is  now  £ dependent) . This  in  turn  is  used  to  solve  the  remaining 
equations  (5),  (6),  (7),  and  (3),  as  well  as  the  shock  relations. 


A significant  aspect  of  the  numerical  integration  is  the  fact 
that  the  shock  wave  is  fixed  as  part  of  the  transformation  from  the 


To  illustrate  these  results  we  consider  a 10%  parabolic  profile  and 

an  upstream  Mach  number  of  five.  The  following  figure  shows  the  map  of  the 

previous  figure  to  the  physical  plane. * equal  constant  are  streamlines 

. ♦ 

(not  all  of  which  are  drawn),  and  Q equal  oonstant  are  C characteristics. 


•1/- 


•*  •c 


A complete  solution  in  graphical  form  is  obtained  on  representing 
the  physical  quantities  on  the  streanlines  and  C+  characteristics. 

As  an  illustration  the  next  two  figures  show  the  variation  entropy 
and  r”(*  ,0).  The  last  figure  is  especially  noteworthy  since  it  gives 
a posteriori  support  for  our  initial  assumption,  vis.,  that  r is  a good 
approxiation  constant  on  streamlines. 


i 


A number  of  additional  cases  have  been  worked  out  and  further  study 
of' these  cases  should  be  made. 

A report  on  the  progress  made  so  far  is  now  in  preparation. 
References 

11]  W.D.  Hayes,  and  R.  Probstein,  Hypersonic  Flow  Theory  (2nd  ed.). 
Academic  Press,  Hew  York,  (1966)  p.  497ff. 

12]  3.3.  Mahony,  J.  Aeronaut.  Sci.  22,  673-680  (1955). 

[3]  R.  Meyer,  Quart.  Appl.  Math.  1£,  433-436  (1957). 
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4.  Aerodynamic  Wakes 

We  have  considered  the  development  of  a wake  which  follows  the 
passage  of  a lifting  body.  As  is  well  known,  a persistent  vortical 
structure  known  as  "horseshoe  vortices"  is  found  in  such  a wake.  Our 
efforts  have  been  directed  at  determining  simplified  equations  which 
govern  this  wake  structure.  We  have,  however,  included  compressibility, 
entropy,  and  dissipative  effects  in  the  calculations. 

As  a starting  point  of  the  analysis,  it  is  assumed  that  the  stream- 
wise  direction,  x,  experiences  slow  variation  as  compared  with  variation 
in  the  cross  (Treffitz)  plane,  (y,  z) . The  method  of  treating  this 
problem  is  based  on  two  scale  perturbation  techniques.  To  summarize 
the  results  we  denote  the  crosswise  velocity  by  u^  ■ (v,w  ) and  the 
longitudinal  velocity  by  u . It  is  then  found  that 


(1)  * Ui.  - 0 

where  V Jl  “ 9 $"7!  so  that  a crossw*se  stream  function  j exists,  i.e. 

(u,  W)  - (-  ^ t . £ f) 

The  stream  function  is  then  found  to  satisfy 

(2)  + <U1  * V i>)£>  - $ 7 1 

where  CO  " is  the  Treffitz  plane  vorticity  and  if  is  a viscous 

coefficient.  The  longitudinal  velocity  u also  satisfies  an  equation  of 
the  same  form. 


wV-.. 
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(3) 

In 

(4) 


+ (ux  * VjJJu  - j 

contrast  to  (2)  and  (3)  the  temperature  is  found  to  be  governed  by 

d . • - 2 

+ (uj.  7x)jT  * S 7 xT 


where *5  is  a heat  conductivity  coefficient. 

As  equations  (2)  and  (3)  indicate,  vorticity  and  entropy  effects 
decouple  and  give  rise  to  different  diffusion  rates.  The  total  vorticity 
to  the  present  order  is  given  by 


(*> 


zhi  £u, 

t>  • ' 2 y] 


and  as  a result  the  twisting  and  stretching  of  vortex  lines  occurs 
within  the  framework  of  this  analysis. 

The  work  just  described  is  independent  of  Mach  number  and  is  thus 
equally  valid  for  supersonic  as  well  as  subsonic  flows. 

In  addition  to  the  above  work  we  have  examined  the  effect  of  a real 
atmosphere  on  the  forms  of  the  governing  equation.  Statification 
effects  and  shearing  winds  have  been  included  in  a further  analysis. 

The  resulting  equations  are  again  in  the  form  of  nonlinear  diffusion 
equations— but  are  now  completely  coupled. 

A full  report  of  this  work  appears  in  the  manuscript: 

L.  Sirovich  and  T.H.  Chong,  On  the  Equations  Governing  a Trailing  Wake 
which  has  been  submitted  for  publication. 
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5.  Two-Dimensional  Unsteady  Plow 

Consider  uniform  (upstream)  two-dimensional  inviscid  flow  past  a 
body  which  is  undergoing  unsteady  motion.  Assume  that  the  motion  of 

« 

the  body  is  "slow"  in  the  sense  that  the  reduced  frequency 


is  small. 6)is  a typical  oscillation  frequency  of  the  body  motion,  L 
is  the  body  length,  and  U is  the  supersonic  uniform  upstream  flow.  From  the 
assumed  weakness  of  the  shock  waves,  we  also  assmme  that  the  Riemann 
invariant  along  the  upstream  going  characteristic  takes  a negligible 
jump  across  the  front  shock  (see  Sec.  2).  Then  from  simple  wave 
transonic  perturbation  theory,  we  obtain  the  governing  equations  for 
the  normalized  velocity  perturbations,  (u,  v) , in  the  form 


y * A fu  )Vi  = Zf?/ ^ r 


where 


Vt  ♦ Fir  * G. 


and  M “ U/c  is  the  upstream  unperturbed  Mach  number,  ft  * M -1, 
o 

V - (£+1)m2. 

The  boundary  condition  on  v is 


v* 


-23- 


where  y * (x,  t)  is  the  body  and 

< r-  rf-  cj  r'w/c'cv) 

(3) 

* « -t-  y/fitv) 

From  shock  analysis,  we  find  that  the  shock  surface  is  governed  by 
an  equation  of  the  form 


where  the  functions  A,  B,  and  C are  explicit,  Eqns.  (3)  and  (4)  together 
with  the  solution  for  v,  (2),  define  a shock  surface  in  the  physical  plane. 

For  purposes  of  illustration,  we  further  assume  that  v is  small 
(this  puts  the  analysis  in  the  supersonic  regime),  and  expand  F and  G 
in  powers  of  v,  then  (4)  yields  (up  to  quadratic  terms  in  v)  a solution 
in  the  form 


<r  = •r't.(t') 

'X.  * i>  -*  'C.C'S') 


where  s is  arc  length,  ((T  ) , 'F  (^  ) , y ('d  ) ) parameterizes  a curve 

o w o o 

on  the  shock  surface,  and  the  functions  a and  b are  explicit.  By 
eliminating  parameters  s and  from  (5),  we  obtain  an  approximate  solution 
to  (4)  in  the  form  y ■ y((T»/t)« 

As  a sample  calculation,  we  have  considered  a pulsating  wing 
defined  by 


. ; 
j y > «XO">')C '* 

I 

for  0<  6(04_x<l,  0 < & < 1.  In  this  case  all  steps  outlined  above 
can  be  carried  out  explicitly  to  yield  a shock  surface  as  well  as 
conditions  on  the  moving  body. 

At  the  present  stage  of  development  we  do  not  feel  that  this  work 
warrants  publication. 


In  a similar  manner  two-dimensional  supersonic  theory  was  considered. 
From  this  analysis  we  found  that  Whitham  theory  (1st  order)  with  only 
minor  modification  may  he  extended  to  second  order.  This  theory  still 
places  characteristics  and  shock  waves  incorrectly  to  this  order.  A 
new  equation  which  has  the  form  (in  two  dimensions) 

*<&,*)%  ♦ pC<p,M)q>„'0 

has  been  derived  and  shown  to  be  valid  uniformly  through  the  second  order. 
This  equation  nay  be  integrated  by  standard  means. 

Although  this  work  may  be  of  analytic  interest,  it  is  superseded 
by  the  work  described  in  Sec.  3.  As  a result  we  have  not  pursued  it 
further  and  do  not  contemplate  a manuscript  on  this  topic. 
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